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Summary 

The Joint Army, Navy, NASA, Air Force (JANNAF) rocket 
engine performance prediction procedure is based on the use 
of various reference computer programs. One of the reference 
programs for nozzle analysis is the Two-Dimensional Kinetics 
(TDK) Program. The purpose of this report is to calibrate the 
JANNAF procedure that has been incorporated into the 
December 1984 version of the TDK program for the high- 
area-ratio rocket engine regime. The calibration was accom- 
plished by modeling the performance of a 1030: 1 rocket nozzle 
tested at NASA Lewis Research Center. A detailed description 
of the test conditions and TDK input parameters is given. 

The results indicate that the computer code predicts delivered 
vacuum specific impulse to within 0. 12 to 1.9 percent of the 
experimental data. Vacuum thrust coefficient predictions were 
within ±1.3 percent of experimental results. Predictions of 
wall static pressure were within approximately ± 5 percent 
of the measured values. An experimental value for inviscid 
thrust was obtained for the nozzle extension between area ratios 
of 427.5 and 1030 by using an integration of the measured 
wall static pressures. Subtracting the measured thrust gain 
produced by the nozzle between area ratios of 427.5 and 1030 
from the inviscid thrust gain yielded experimental drag 
decrements of 10.85 and 27.00 N (2.44 and 6.07 lb) for 
mixture ratios of 3.04 and 4.29, respectively. These values 
correspond to 0.45 and 1.11 percent of the total vacuum thrust. 
At a mixture ratio of 4.29, the TDK predicted drag decrement 
was 16.59 N (3.73 lb), or 0.71 percent of the predicted total 
vacuum thrust. 

Introduction 

In 1975, the Joint Army, Navy, NASA, Air Force 
(JANNAF) Rocket Engine Performance Working Group 
developed and documented a methodology to model rocket 
engine systems. This methodology, outlined in reference 1, 
was developed to create an industry and government reference 
for rocket engine performance prediction. 

The JANNAF prediction procedure makes use of various 
reference computer programs. One of the reference programs 
is the Two-Dimensional Kinetics (TDK) Program for nozzle 
analysis (ref. 2). TDK was originally developed under the 
auspices of the JANNAF working group. At that time, TDK 


performed a two-dimensional, inviscid calculation of rocket 
nozzle performance. Over the years, the TDK program has 
been extended to include a prediction of viscous effects on 
nozzle performance using the JANNAF procedure. 

When the JANNAF procedure was developed, large-area- 
ratio rocket nozzles extended to area ratios of 100. With the 
recent effort to develop engines for applications such as the 
orbital transfer vehicle, rocket nozzle designs with area ratios 
of 1000 or larger are being examined. Because these high- 
area-ratio nozzles create a new performance prediction 
domain, it is unclear how well the JANNAF procedure will 
predict. Therefore, there is a need to calibrate the procedure 
for this rocket engine regime. 

The purpose of this report is to calibrate the JANNAF 
procedure that was incorporated into the December 1984 
version of the TDK program. The calibration is accomplished 
by modeling the performance of a 1030: 1 rocket nozzle tested 
at NASA Lewis Research Center. A detailed description of 
the test conditions and TDK input parameters is given. 

This report presents experimental vacuum thrust and vacuum 
specific impulse I spV data for an optimally contoured nozzle, 
which was extended to an exit area ratio of 1030 and could 
be truncated to an exit area ratio of 427.5. The nozzle was 
tested using a gaseous hydrogen and gaseous oxygen 
combustion system at a nominal chamber pressure of 2413 
kN/m 2 (350 psia) and a propellant mixture ratio OIF range 
of 2.78 to 5.49 (ref. 3). The experimental thrust and l spV 
results are compared to the theoretical predictions obtained 
from the TDK computer code. 

Experimental wall static pressures were used to quantify the 
inviscid thrust gain between the area ratios of 427.5 and 1030. 
By comparing this inviscid thrust to the measured thrust gain, 
we obtained a value for the shear (or drag) force. Corre- 
sponding values were obtained from the TDK program and 
compared to the experimental results. 


Background 

Test Facility 

Testing was done in the new altitude test capsule at the 
NASA Lewis Rocket Engine Test Facility (RETF). Figure 1 
is an illustration of RETF with cutaway views of the test 
capsule and spray cooler. The operation of the facility was 
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Figure 1.— Rocket Engine Test Facility (RETF) with cutaway views of test capsule and spray cooler. 


as follows. When the engine was fired, the exhaust gases 
flowed into the diffuser where the kinetic energy of the exhaust 
was used to accomplish some of the altitude pumping. From 
the diffuser, the exhaust gases flowed into the spray cooler, 
where approximately one-half of the exhaust was condensed 
to water and flowed out the vertical drain line. The other half 
was pumped by the gaseous nitrogen ejectors shown mounted 
on top of the spray cooler. The pressure obtained in the test 
capsule was in the range of 0.207 to 0.34 kN/m 2 (0.03 to 


0.05 psia). A more indepth description of the facility and test 
apparatus can be found in reference 3. 

The thrust stand used in this facility was capable of 
measuring 13.34 kN (3000 lb) full scale and was attached to 
a foundation that was separate from the test capsule bulkhead. 
The thrust stand was designed to have a 2<r (standard deviation) 
variation of less than ±0.1 percent of full scale, and it was 
calibrated against a load cell that had a 2<r variation of less 
than ±0.05 percent of full scale. 



Figure 2.— Altitude test capsule— 1030: 1 nozzle being installed in thrust stand. 
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Test Hardware 

Figure 2, which is a photograph of the inside of the test 
capsule, shows the thrust stand with the 1030:1 nozzle in the 
process of being installed. The injector used during these tests 
had 36 gaseous oxygen shower-head elements, with the 
gaseous hydrogen flowing through a porous face plate. The 
solid copper combustion chamber was 15.24 cm (6 in.) long 
and had an inside diameter of 5.2197 cm (2.055 in.). It was 
uncooled (heat-sink) hardware that relied on its thermal 
capacitance to survive the short firings (<3 sec). 

The 1030: 1 nozzle was made of three sections. The first, 
the throat section, connected to the combustion chamber, 
converged to the 2.54 cm (1 in.) throat, and then diverged 
to an area ratio of 29.9. The throat section was made of nickel 
with a ceramic coating on the inner wall and had a water jacket 
around the throat that extended to an area ratio of 
approximately 5. The next section of the nozzle expanded to 
an area ratio of 427.5 and was made of 0.635 cm (0.25 in.) 
carbon steel. The last section of the nozzle extended to an area 
ratio of 1030 and was also made of 0.635-cm (0.25-in.) carbon 
steel. The carbon steel sections were heat-sink hardware and 
had no ceramic coatings. 

The 1030:1 nozzle contour was designed using the Rao 
Nozzle Contour Program (ref. 4) and the Boundary-Layer 
Integral Matrix Procedure (BLIMP- J) (ref. 5) program. The 
Rao program provided the inviscid optimal nozzle contour, 
and the BLIMP-J program estimated the boundary-layer 
displacement thickness along the nozzle. The combination of 
the two results provided the nozzle coordinates that were used 
to make the hardware. Figure 3 shows the nozzle contour and 
coordinates. A further discussion of the nozzle design is 
presented in reference 3. 

Analysis 

The TDK computer program evaluates the two-dimensional 
nonequilibrium chemistry and viscous effects on the 
performance of rocket exhaust nozzles. Version 2.4 
(December 1984) of TDK was used in this report. TDK 
consists of a master control module (MCM) and five compu- 
tational modules: ODE, ODK, TRAN, MOC, and BLM. The 
MCM controls the execution of TDK and processes the output. 

The ODE (one-dimensional equilibrium) module calculates 
ideal engine performance assuming either chemical equilibrium 
composition or a frozen chemical composition at rocket 
chamber conditions. ODE uses the free-energy minimization 
method to compute equilibrium conditions for any assigned 
enthalpy and pressure. 

The ODK (one-dimensional kinetics) module calculates the 
inviscid, one-dimensional nonequilibrium nozzle expansion of 
gaseous propellant exhaust. 

The TRAN (transonic flow) module uses the chemical 
information computed by ODK to estimate two-dimensional 
effects in the transonic region of the nozzle throat. The purpose 
of these calculations is to approximate an initial data line across 


the nozzle throat in order to start the method of characteristics 
(MOC) calculations. 

The MOC module uses the method of characteristics to 
construct a finite-difference mesh by tracing gas streamlines 
and characteristic surfaces. By determining the properties of 
the exhaust at the mesh points, MOC is able to calculate the 
loss in nozzle performance caused by flow divergence. 

The BLM (boundary-layer module) calculates compressible 
laminar and turbulent wall boundary layers in axisymmetric 
nozzles. BLM uses the two-point finite-difference method 
developed by Keller and Cebeci (ref. 6) to calculate the 
boundary-layer properties and Cebeci-Smith eddy-viscosity 
formulation (ref. 7) to model the turbulent boundary layer. 

The experimental hardware specifications and test conditions 
were used to write the input files to TDK so that this program 
could accurately model the nozzle performance. The input 
variables that described the nozzle inlet geometry are listed 
in table I. The nozzle contour coordinates that were used are 
shown in figure 3. Table II shows the experimental results that 
were also used in the TDK input files: effective chamber 


NOZZLE COORDINATES 


AXIAL DISTANCE 
FROM THROAT 

RADIUS J 

cm 

in. 

cm 

in. 

0.0000 

0.0000 

1.2700 

0.5000 

.3929 

.1547 

1.4371 

.5658 

.4641 

.1827 

1 . 4961 

.5890 

.6068 

.2389 

1.6190 

.6374 

.7503 

.2954 

1.7404 

.6852 

.8230 

.3240 

1.8031 

.7099 

1.3246 

.5215 

2.2426 

.8829 

1.7844 

.7025 

2.6515 

1.0438 

2.3777 

.9361 

3.1643 

1.2458 

3.2062 

1.2623 

4.2001 

1.6536 

7.0256 

2.7660 

6.6703 

2.6261 

7.8931 

3.1075 

7.2426 

2.8514 

9.6269 

3.7901 

8.3320 

3.2803 

10.6505 

4.1931 

8.9433 

3.5210 

11.6738 

4.5%0 

9.5341 

3.7536 

12.9022 

5.07% 

10.2189 

4.0232 

15.3429 

6.0405 

11.5108 

4.5318 

16.5392 

6.5115 

12.1150 

4.7697 

19.5651 

7.7028 

13.5702 

5.3426 

23.3688 

9.2003 

15.2710 

6.0122 

25.4869 

10.0342 

16.1651 

6.3642 

29.5410 

11.6303 

17.7871 

7.0028 

33.7297 

13.2794 

19.3558 

7.6204 

36.29% 

14.2912 

20.2705 

7.9805 

38.86% 

15.3030 

21.1524 

8.3277 

41.4193 

16.3068 

21.9977 

8.6605 

47.2194 

18.5903 

23.8201 

9.3780 

51.1703 

20.1458 

24.9895 

9.8384 

55.1213 

21.7013 

26.1064 

10.2781 

60.4944 

23.8167 

27.5486 

10.8459 

71.1091 

27.9957 

30.1694 

11.8777 

76.2211 

30.0083 

31.3365 

12.3372 

90.63% 

35.6849 

34.3444 

13.5214 

105.0371 

41.3532 

36.9933 

14.5643 

113.0838 

44.5212 

38.3365 

15.0931 

128.5725 

50.6191 

40.6598 

16.0078 
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Figure 3.— Nozzle contour and coordinates. 
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TABLE I.— TWO-DIMENSIONAL KINETICS (TDK) INPUT 
VARIABLES 


Parameter 

TDK variable 

Value 

Throat radius, cm (in.) 

RSI 

1.27(0.5) 

Inlet contraction ratio 

ECRAT 

4.223 

Inlet wall radius 3 

RI 

2.0 

Inlet angle, deg 

THETAI 

25.0 

Upstream wall radius 
of curvature 3 

RWTU 

2.0 

Downstream wall radius 
of curvature 3 

RWTD 

0.4 

Nozzle attachment angle, deg 

THETA 

39.41 

Nozzle exit angle, deg 

THE 

b 7.94 


formalized by throat radius. 
t’THE = 15.5 for truncated contour. 


TABLE II.— EXPERIMENTAL RESULTS 


Reading 

Nozzle exit 
expansion 
area ratio, 
€ 

Effective 
chamber pressure 
Pc.e 

Propellant 

mixture 

ratio, 

O/F 

Fuel injection 
pressure 

Fuel injection 
temperature 

Oxidizer 

injection 

pressure 

Oxidizer 

injection 

temperature 

Propellant flow 
rate 

kN/m 2 

psia 

K 

•R 

kg/sec 

lb/sec 

kN/m 2 

psia 

kN/m 2 

psia 

K 

“R 

112 


2482 


3.84 

3061 


285.6 

514.1 


407.4 

279.2 

502.5 

0.5266 

1.161 

113 



2461 

356.9 

4.36 

2956 


284.2 

511.5 

2818 

408.7 

277.0 

498.6 

.5334 

1.176 

114 



2488 


5.08 

2912 

Q 

283.9 

511.1 

2890 

419.2 

275.8 

496.4 

.5543 

1.222 

115 




355.3 

5.49 

2843 

412.4 

283.5 


2870 

416.3 

275.9 

496.6 

.5552 

1.224 

117 



2456 

356.2 

3.19 

3160 

458.3 

281.1 

506.0 

2735 

396.7 

275.1 

495.2 

.5094 

1.123 




2449 

355.2 


2950 

427.8 

294.4 



406.6 

287.5 

517.5 

.5316 

1.172 

121 



2482 


4.11 

3013 

437.0 


531.0 

2832 

410.7 

288.3 

518.9 

.5293 

1.167 

123 



2449 

355.2 

3.19 

3152 

457.2 

295.6 

532.0 

2732 

396.2 

289.4 

521.0 

.5039 

1.111 

124 



2492 

361.4 

2.78 

3328 

482.7 

295.9 

532.6 

2752 

399.1 

289.5 

521.1 

.5058 

1.115 

125 



2441 

354.0 

3.74 

3028 

439.1 

296.3 


2763 

400.7 

289.2 

520.6 

.5135 

1.132 

136 

427.5 

2383 

345.6 


3123 

452.9 

292.1 

525.7 

2648 

384.1 

287.6 

517.7 

.4917 

1.084 

137 

427.5 

2457 

356.8 

4.29 

2965 

EFIwl 

291.2 

524.2 

2813 

wm 

285.0 

531.0 

.5330 

1.175 


Reading 

Vacuum thrust, 
Fy 

Ambient pressure 
around nozzle, 
Pa 

Characteristic 
exhaust velocity, 
C* 

Characteristic 
exhaust velocity 
efficiency, 

VC"’ 

percent 

Measured 
vacuum thrust 
coefficient, 
c F,y 

Vacuum 

thrust coefficient 
efficiency, 
VCpy 
percent 

Vacuum specific 
impulse, 

hp,V< 

sec 

Vacuum specific 
impulse efficiency, 

1,1 spy 
percent 

N 

lb 

kN/m 2 

psia 

m/sec 

ft/sec 

wm 

2422 

544.4 



2424 

7953 

96.4 

mm 

97.3 

468.9 

92.9 

mm 

2409 

541.6 

.2592 

IB 

2366 

7762 

95.8 


95.3 

460.4 

90.5 

114 

2457 

552.3 

.2530 

mm 

2305 

7562 

95.3 


94.0 

451.9 

88.8 

115 

2448 

550.4 

.2530 

WM 

2261 

7418 

94.9 


93.7 

449.7 

88.3 

117 

2364 

531.5 

.2461 

.0357 

2473 

8115 

97.2 


98.4 

473.4 

94.9 

120 

2429 

546.1 

.2544 

.0369 

2382 

7815 

95.9 


96.0 

466.1 

91.8 

121 

2459 

552.9 

.2654 

.0385 

2403 

7885 

96.2 


96.6 

473.6 

93.5 

123 

2377 

534.3 

.2592 

.0376 

2477 

8128 

97.2 

1.881 

97.7 

481.1 

96.2 

124 

2336 

536.4 

.2441 

.0354 

2502 

8208 

97.8 

1.857 

97.8 

481.3 

97.3 

125 

2406 

541.0 

.2420 

.0351 

2433 

7984 

96.5 

1.912 

97.4 

477.8 

94.6 

136 

2228 

500.9 

.3916 

.0568 

2487 

8158 

97.4 

1.823 

96.0 

462.3 

93.6 

137 

2365 

531.6 

.3916 

.0568 

2384 

7822 

96.0 

1.877 

94.8 

452.6 

90.2 
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pressure, propellant mixture ratio, fuel injection temperature, 
and oxidizer injection temperature. The experimentally 
determined wall temperatures were used and are listed in the 
input files in appendix A. 

Boundary-layer edge conditions and wall temperatures 
within the combustion chamber and convergent nozzle had to 
be estimated because no experimental data were available. 
These values appear in the BLM namelist in the input files. 

The TDK program requires the location of the transition 
from laminar to turbulent flow as input. Therefore, a study 
was performed to determine the approximate boundary-layer 
transition point within the nozzle. By comparing the experi- 
mental heat flux data to predicted heat flux, we determined 
that the boundary layer was laminar over the entire nozzle 
(ref. 8). To model the laminar flow, we instructed the program 
to place the transition point beyond the exit plane of the nozzle. 
For comparison, cases were also run for a boundary layer that 
transitioned to turbulent within the combustion chamber. 

Appendix A contains the TDK input files for experimental 
readings 112 to 115, 120, 121, and 137. Only these experi- 
mental readings could be modeled because the TDK program 
could not run to completion for mixture ratios below 3.84. 
This version of the TDK program was originally written for 
much lower area ratio rocket nozzle conditions. The program 
was unable to run below an O/F of 3.84 because of the low- 
pressure/low-temperature conditions that are predicted as the 
flow expands to an area ratio of 1030. 

The program was instructed to calculate the boundary-layer 
displacement thickness for the actual nozzle contour and to 
use it to obtain the displaced or inviscid contour. This inviscid 
contour was then run through the MOC module to obtain the 
final predictions of I spV and thrust. 

Discussion of Results 

This section presents the analytical and experimental results 
for an optimally contoured nozzle which expanded to an area 
ratio of 1030 and truncated to an exit area ratio of 427.5. The 
analytical results predicted by the December 1984 version of 
TDK for the experimental readings 1 12 to 1 15, 120, 121 , and 
137 are presented in table HI. The corresponding experimental 
results are presented in table II. 

By evaluating the measured heat flux from the nozzle and 
the estimated flow conditions within the nozzle, we determined 
that the boundary layer behaved as a laminar boundary layer 
throughout the entire nozzle (ref. 8). Thus, the TDK program 
was instructed to assume laminar flow in determining viscous 
effects. To learn the effect on predicted performance, we also 
made separate runs of TDK for a boundary layer that 
transitioned to turbulent within the combustion chamber. 

In order to predict the delivered l sp y, TDK predictions 
must be adjusted to account for energy-release losses. Energy - 
release losses consist of two parts: vaporization losses and 
mixture ratio distribution losses. Vaporization losses are due 


TABLE m — TDK RESULTS 


Reading 

Nozzle exit 
expansion 
area ratio, 
€ 

Effective 

chamber 

pressure, 

P ce 

Measured 

propellant 

mixture 

ratio, 

O/F 

Predicted 
propellant 
flow rate 

kg/sec 

Ib/sec 

kN/m 2 

psia 

112 

1030 

2482 

360.0 

3.84 

0.5034 

1.1099 

113 



2461 

356.9 

4.36 

.5065 

1.11667 

114 



2488 

360.9 

5.08 

.5245 

1.15628 

115 



2450 

355.3 

5.49 

.5243 

1.15595 

120 



2449 

355.2 

4.30 

.5027 

1.10834 

121 



2482 

360.0 

4.11 

.5066 

1.1169 

137 

427.5 

2457 

356.8 

4.29 

.5051 

1.11352 


Reading 

Computer code 

TDK/BLM, laminar 

Predicted 
characteristic 
exhaust velocity, 

c- 

Predicted vacuum 
thrust, 

F, 

Predicted 

vacuum 

thrust 

coefficient, 

r 

'-'F.V 

Predicted 

thrust 

coefficient 

efficiency, 

1C F y 

percent 

N 

lb 

m/sec 

ft/ sec 

112 

2502.37 

8209.89 

2383.18 

535.76 

1.8917 

96.09 

113 

2466.65 

8092.70 

2393.78 

538.144 

1.916 

95.37 

114 

2408.94 

7903.35 

2467.18 

554.644 

1.9527 

94.56 

115 

2373.03 

7785.52 

2452.14 

551.263 

1.9708 

93.98 

120 

2473.29 

8114.47 

2378.05 

534.607 

1.9125 

95.42 

121 

2487.25 

8160.27 

2399.67 

539.468 

1.9044 

95.75 

137 

2473.41 

8114.85 

2332.30 

524.321 

1.8669 

94.23 


Reading 




Computer code 



ODE 

ODK 

MOC 


TDK/BLM, laminar 


Predicted vacuum specific impulse, 
hp, V’ 
sec 

Predicted vacuum 
specific impulse 
(adjusted), 
hp.V' 

sec 

Predicted vacuum, 
specific impulse 
efficiency (adjusted), 

percent 

112 

504.8 

498.45 

494.211 

482.714 

465.34 

92.18 

113 

507.4 

498.44 

494.622 

481.921 

461.68 

90.99 

114 

509.2 

496.89 

493.59 

479.678 

457.13 

89.77 

115 

509.5 

494.77 

491.734 

476.894 

452.57 

88.83 

120 

507.8 

499.02 

495.126 

482.347 

462.57 

91.09 

121 

507.0 

499.21 

495.287 

483.006 

464.65 

91.65 

' 137 

502.0 

494.19 

480.215 

470.87 

452.035 

90.05 


to incomplete liquid droplet vaporization at the nozzle throat. 
Mixture ratio distribution losses are due to nonuniform 
distribution of the vaporized propellant at the nozzle throat. 
As described in reference 9, experimental characteristic 
exhaust velocity efficiency t \ c • can be used as an estimate of 
specific impulse energy release losses. Thus, the TDK 
predictions were multiplied by the experimental characteristic 
velocity efficiencies shown in table II to account for energy- 
release losses. These results are labeled “Adjusted 
TDK/BLM” on figures 4, 5, and 8. 

Performance Results— 1030:1 Area Ratio Nozzle 

The basic measure of rocket engine performance is specific 
impulse I sp . In figure 4, the predicted thrust chamber losses 
from ideal or maximum performance are presented. The 
ODE curve represents the predicted ideal, one-dimensional 
equilibrium values of I sp y. The ODK curve indicates the 
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predicted results for one-dimensional, nonequilibrium flow. 
Thus, the drop in I sp V from ODE to ODK represents the loss 
in performance due to kinetics. For the 1030:1 nozzle, these 
losses are estimated to be 1.3 to 2.89 percent of maximum 
l sp y over the mixture ratio range from 3.84 to 5.49. 

Points on the MOC curve are obtained from the MOC 
module and represent the inviscid, two-dimensional, nonequili- 
brium predictions. The difference between the ODK and MOC 
curves is the loss in performance due to nozzle divergence 
shape and exit angle. (The actual nozzle contour was used in 
these MOC calculations.) The estimated divergence losses 
range from 0.60 to 0.84 percent. 

The TDK/BLM curve was obtained using the final results 
from the TDK program. These results contain the predicted 
boundary-layer losses from the BLM and MOC calculations 
for the displaced, or inviscid, contour. The difference between 
the MOC curve and the TDK/BLM curve is the performance 
loss due to laminar boundary-layer effects and is estimated 
to be 2.3 to 2.9 percent maximum I spV for the specified 
mixture ratio range. 

As mentioned previously, the TDK values of I spV were 
adjusted to account for energy release losses. These values 
are shown on the adjusted TDK/BLM curves. From a point- 
to-point comparison, the adjusted TDK predictions for a 
completely laminar boundary layer are within 0.3 to 1.9 
percent of the experimental readings modeled. Based on the 
results shown in figure 4, there appears to be no correlation 
between the accuracy of the prediction and the mixture ratio 
at which the prediction is made. The adjusted predictions for 
a turbulent boundary layer are 2.3 to 5.0 percent lower than 
the experimental results. Overall, the turbulent predictions are 
approximately 2.5 percent lower than the laminar predictions. 
This amounts to roughly a 15-sec drop in I spV . Thus, proper 
determination of the boundary-layer characteristics prior to 
making a performance prediction can be important. 

Figure 5 is a plot of the thrust chamber performance 
efficiency. Efficiency is calculated by dividing the 
experimental and TDK results by the ideal (ODE) values. As 
shown, the performance efficiency increases as mixture ratio 
decreases. The experimental and predicted values of efficiency 
compare to the same degree as in figure 4. 

The vacuum thrust coefficient C F V is a quantity that 
reflects the design quality of a nozzle. It is an indication of 
the thrust augmented by the gas expansion through the nozzle 
as compared with the thrust that would be generated if the 
chamber pressure acted over the throat area only. In figure 6, 
the experimentally obtained values of C FV are presented 
along with the TDK predictions for a completely laminar 
boundary layer and the TDK predictions for a turbulent 
boundary layer. As indicated in the plot, thrust produced by 
the nozzle increased as mixture ratio increased. The difference 
between the experimental and TDK/BLM laminar results is 
within ±1.3 percent. For a turbulent boundary layer, the 
predictions fall approximately 3.5 percent below the 
experimental results. 




PROPELLANT MIXTURE RATIO. 0/F 

Figure 5.— Thrust chamber performance efficiency. Area ratio, e, 1030. 
Specific impulse efficiency i} ; is based on ideal one-dimensional 
equilibrium (ODE) results. 


O EXPERIMENTAL RESULTS 



1.80 1 1 ' 1 1 1 1 1 

2.5 3.0 3.5 A.O A. 5 5.0 5.5 6.0 

PROPELLANT MIXTURE RATIO, 0/F 

Figure 6.— Experimental and predicted vacuum thrust coefficient. Area ratio, 
6, 1030. 
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Figure 7 is a plot of the vacuum thrust coefficient efficiency 
based on ideal results for the 1030:1 nozzle. As shown 
previously with performance efficiency, vacuum thrust effi- 
ciency increases as the mixture ratio decreases. These curves 
compare in the same manner as those in figure 6. 

Performance Results— Truncated Contour 

The TDK program was used to model the experimental 
results for the truncated contour and the test conditions in 
reading 137. For I spV the program predicted 1.6-percent 
kinetic losses, 2.8-percent divergence losses, and 1.9-percent 
laminar boundary-layer losses. After adjusting the TDK/BLM 
laminar l sp ,v prediction for energy release losses, the 
prediction of 452.04 sec was within 0.12 percent of the 
experimental value of 452.6 sec. A turbulent boundary-layer 
assumption for this configuration yields an I sp V of 441.28 sec 
(2.5 percent lower than the experimental I sp y). 

The predicted C FV for the truncated contour with a 
completely laminar boundary layer is 1 .8669 and is within 0.5 
percent of the experimental value of 1.877. The TDK 
prediction for the truncated contour with a turbulent boundary 
layer yields a C F V that is 2.9 percent lower than the 
experimental C F V . 

The nozzle was truncated to experimentally determine the 
performance gain between the area ratios of 427.5 and 1030 
and to validate predictions of performance over the length of 
the nozzle. Figure 8 is a plot of the predicted I sp y over the 
length of the 1030:1 nozzle for a mixture ratio of 4.3. The 
two points plotted indicate the experimental values of I spV 
obtained from the full and truncated contours. As shown, there 
is very good agreement between the TDK laminar predictions 
and the experimental results. The TDK turbulent predictions 
fall approximately 10 to 20 sec below the laminar curve 
between an area ratio of 100 and the exit. 

Pressure Integration Results 

Figure 9 is a plot of the wall static pressure distribution for 
the 1030: 1 nozzle. It is expressed as the ratio of wall static 
pressure P w to effective chamber pressure P c e for reading 
115. Effective chamber pressure is an estimate of nozzle throat 
total pressure. The method used to obtain P c e is discussed in 
reference 3. There are two significant observations to be drawn 
from this figure. The first observation indicates whether the 
nozzle is flowing full or whether the flow has separated from 
the wall. For the experimental data reported, there was no 
separation, as the pressure distribution continued to expand 
all the way to the exit plane of the nozzle. The second 
observation compares experimental and analytically predicted 
values of pressure. TDK predicted static pressures are within 
±5 percent over most of the nozzle. The greatest difference 
between prediction and experiment is at the area ratios of 12 
and 1000. At these points, the experimental values are 
approximately 15 percent higher. This could be an indication 



Figure 7.— Efficiency of thrust amplification based on one-dimensional 
equilibrium (ODE) results. Area ratio, e, 1030. 



Figure 8.— Predicted and experimental thrust chamber performance over length 
of nozzle. Area ratio, e, 1030. 



Figure 9.— Static wall pressure distribution for extended nozzle. Reading 115; 
area ratio, e, 1030; propellant mixture ratio, O/F, 5.49; effective chamber 
pressure, P ce , 2450 kN/m 2 (355.3 psia). 
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AREA RATIO. £ 

Figure 10.— Static wall pressure distribution for truncated nozzle. Reading 
137; area ratio, e, 427.5; propellant mixture ratio, O/F, 4.29; effective 
chamber pressure, P ce , 2460 kN/m 2 (356.8 psia). 

of separated flow near the throat and of boundary-layer 
feedback at the exit. Both of these phenomena would result 
in a higher pressure measurement than predicted. Further 
investigation into this matter is beyond the scope of this report. 

Figure 10 is the wall static pressure distribution for the 
truncated nozzle (reading 137). The TDK predicted pressures 
are within ±5 percent of experimental pressures, except for 
area ratios of approximately 300 to the exit. Over this part 
of the nozzle, the experimental values are approximately 15 
percent higher than predicted values. Because the result of this 
comparison was not the same as that for the full contour, a 
comparison was made of the pressures measured in readings 
113 (full contour) and 137 (truncated contour). In both cases, 
P ce and O/F were similar. The comparison indicated that the 
pressures measured at the same area ratios agreed to within 
±0.5 percent up to an area ratio of 200. For the pressures 
measured at area ratios of 200, 300, and 388, the values were 
5 to 10 percent higher for reading 137. This indicates that 
truncating the nozzle affected the experimental static pressure 
measurements, which may explain the 15-percent difference 
between experimental and theoretical values in figure 10. The 
higher pressures could be occurring as a result of the response 
of the subsonic boundary layer to the change in capsule 
pressure from reading 113 to reading 137. Further testing is 
needed to study this phenomenon. 

A calculation was performed using the measured wall static 
pressure distribution for the 1030:1 nozzle. The area under 
the pressure plot was determined from an area ratio of 427.5 
to 1030. This pressure integration yields the experimental 
inviscid thrust gain achieved with the addition of the nozzle 
extension from 427.5:1 to 1030:1. When the actual thrust gain 
(the difference in measured thrust for the full and truncated 
contours) is subtracted from the inviscid thrust gain, the result 


is the shear or drag force that is produced between the area 
ratios of 427.5 and 1030. A detailed discussion of the pressure 
integration procedure appears in appendix B. 

Integration of the measured static pressures between the area 
ratios of 427.5 and 1030 yielded values of 80.07 N (18.0 lb), 
or 3.35 percent of the total vacuum thrust, at a mixture ratio 
of 3.04 and 82.69 N (18.59 lb), or 3.40 percent of the total 
vacuum thrust, at a mixture ratio of 4.29. The measured thrust 
gain at those mixture ratios was 69.21 N (15.56 lb), or 2.90 
percent, and 55.69 N (12.52 lb), or 2.29 percent, respectively. 
Thus, the shear or drag force was as follows: 



Mixture ratio 

3.04 

4.29 

Inviscid thrust gain, N (lb) 

80.07 (18.00) 

82.69 (18.59) 

Actual thrust gain, N (lb) 

69.21 (15.56) 

55.69 (12.52) 

Drag, N (lb) 

“10.86 (2.44) 

“27.00 (6.07) 

Vacuum thrust, percent of total 

0.45 

1.11 


a Actual thrust gain subtracted from inviscid thrust gain. 

The corresponding TDK predictions at a mixture ratio of 4.29 
were 62.34 N (14.015 lb) inviscid thrust, or 2.67 percent of 
the total vacuum thrust, and 45.75 N (10.286 lb), or 1.96 
percent of the thrust gain. Therefore, the predicted shear or 
drag force was as follows: 


Mixture ratio 

4.29 

Inviscid thrust gain, N (lb) 

62.34 (14.02) 

Actual thrust gain, N (lb) 

45.75 (10.29) 

Drag, N (lb) 

“16.59 (3.73) 

Vacuum thrust, percent of total 

0.71 


a Actual thrust gain substracted from inviscid thrust gain. 

Thus, the TDK prediction of drag was lower than the 
experimental value by 10.23 N (2.3 lb), or 0.4 percent of the 
total vacuum thrust. 

Figure 1 1 was developed by extending the examination of 
inviscid thrust and drag decrement. In this figure, the predicted 



AREA RATIO, £ 

Figure 11.— Predicted thrust gain over length of nozzle. Reading 115; 
area ratio, e, 1030; propellant mixture ratio. O/F, 5.49. Nondimensional 
axial distance is axial distance divided by throat radius. 
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drag force and predicted inviscid thrust produced per unit 
length of axial distance are plotted for the 1030: 1 nozzle at 
a mixture ratio of 5.49 (reading 115). The distance between 
the curves represent the actual thrust gain per unit length for 
a specific area ratio. An analysis of this kind could be used 
to examine the advantage of either extending a nozzle contour 
or truncating it. As the need for higher area ratio nozzles 
develops, the examination of thrust gain with area ratio will 
become more important. 

Summary of Results 

Predictions from the December 1984 version of the TDK 
nozzle analysis program were compared to experimental 
results of a test-fired rocket nozzle. The hardware tested was 
a heat-sink nozzle optimally expanded to an area ratio of 1030 
and designed so that it could be truncated at an area ratio of 
427.5. Test conditions, which included mixture ratio, propel- 
lant (gaseous oxygen and gaseous hydrogen) temperatures, 
chamber pressure, and nozzle wall temperatures, were used 
in the computer program to accurately model the nozzle 
performance. 

An evaluation of the measured heat flux from the nozzle 
and the estimated flow conditions within the nozzle determined 
that the boundary layer behaved as a laminar boundary layer 
throughout the entire nozzle. Thus, the TDK program was 
instructed to assume laminar flow in determining viscous 
effects. To learn the effect on predicted performance, we also 
made separate runs of TDK for a boundary layer that 
transitioned to turbulent within the combustion chamber. The 
TDK predictions indicated that a 2.5-percent difference in 
vacuum specific impulse can result from using a turbulent 
boundary-layer assumption, instead of a completely laminar 
boundary-layer assumption. Thus, proper determination of the 
boundary-layer characteristics prior to making a performance 
prediction can be important. 


The TDK predictions were compared to the experimental 
results for the full and truncated contours. The parameters that 
were compared were vacuum specific impulse, vacuum 
specific impulse efficiency, vacuum thrust coefficient, vacuum 
thrust coefficient efficiency, wall static pressure, and thrust 
gain between the area ratios of 427.5 and 1030. The results 
of the comparison for the mixture ratio range from 3.84 to 
5.49 were as follows: 

1. The TDK predictions of delivered vacuum specific 
impulse were within 0.12 to 1.9 percent of the experimental 
results. There appears to be no correlation between the 
accuracy of the prediction and the mixture ratio at which the 
prediction was made. The experimental and predicted values 
of efficiency compare to the same degree. 

2. Vacuum thrust coefficient predictions were within ±1.3 
percent of experimental results. Again, the predictions of thrust 
coefficient efficiency compare to the same degree. 

3. The predictions of wall static pressure were within ±5 
percent of experimental results, except at the lowest and 
highest area ratios at which they were measured. At these area 
ratios of 12 and — 1000, the predictions were 15 percent lower 
than measured results. 

4. An experimental value for inviscid thrust was obtained 
for the nozzle extension between area ratios of 427.5 and 1030 
by using an integration of the measured wall static pressures. 
Subtracting the measured thrust gain produced by the nozzle 
between area ratios of 427.5 and 1030 from the inviscid thrust 
gain yielded experimental drag decrements of 10.86 and 
27.00 N (2.44 and 6.07 lb) for mixture ratios of 3.04 and 4.29, 
respectively. These values correspond to 0.45 and 1.11 percent 
of the total vacuum thrust. At a mixture ratio of 4.29, the TDK 
predicted drag decrement was 16.59 N (3.73 lb), or 0.71 
percent total vacuum thrust. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 6, 1987 


Appendix A 
TDK Input Files 


Input for Reading 112 



LOW T 
H 

100 . 

CPHS 

4 . 968 

2 

1 

200 . 

4.968 


2 

H2 
100 . 

6.729 

2 

1 

200 . 

6.560 


2 

H20 
100 . 

7.961 

2 

1 

200. 

7.969 


2 

0 

100 . 

5.665 

2 

1 


9 


5.433 


100. 7.798 1 

200. 7.356 2 

02 2 

100. 6.956 1 

200. 6.961 2 

END LOW T CPHS 

TITLE HIGH E NOZZLE STUDY: E=1000 AND READING = 112 
DATA 
EDATA 

ODE= 1 , ODK=l , TDK= 1 , BLM= 1 , IRPEAT=2 , IOFF=6 , 

RSI=0 . 5 , ASUB=3 . , 1 . 5 , NASUB=2 , IRSTRT=0 , 

ASUP= 1.5,2.0,30.0,200.0,400.0,600.0,1024.0, NASUP=7 , 
ECRAT=4 . 223 , RI=2 . ,THETAI=25. ,RWTU=2. , 

ITYPE=0 , IWALL=4 ,RWTD=0 . 4 ,THETA=39 . 41 , 

THE=7 . 94 , NWS=36 , 

RS= 0.0000, 1.1316, 1.1780, 1.2748, 

1.3704, 1.4198, 1.7658, 2.0876, 


RS= 0.0000, 1.13 

1.3704, 1.4198 

2.4916, 3.0372, 

6.5606, 7.0420 

9.0636, 9.5394, 

12.7284, 14.0056, 

16.6554, 17.3210, 

21.6918, 23.7554, 

29.1286, 30.1862, 3 

ZS= 0.0000, 0.309 

0.5908, 0.6480 

1.8722, 2.5246, 

7.5802, 8.3862 

12.0810, 13.0230, 

20.0684, 23.2606, 

30.6060, 32.6136, 

47.6334, 55.9914, 

82.7063, 89.0425, 

EEND 

REACTANTS 

H 2. 

0 2 . 


)0, 1.1316, 1.1780, 1.: 

t, 1.4198, 1.7658, 2.08; 

3.0372, 5.2522, 5.701 

>, 7.0420, 7.5072, 8.041 

9.5394, 10.6852, 12.024- 

14.0056, 15.2408, 15.9611 

17.3210, 18.7560, 19.6768, 

23.7554, 24.6744, 27.0428 

30.1862, 32.0156, 

), 0.3094, 0.3654, 0.4' 

i, 0.6480, 1.0430, 1.40! 

2.5246, 5.5320, 6.21! 

>, 8.3862, 9.1920, 10.15 

13.0230, 15.4056, 18.4001 

23.2606, 26.5588, 28.582- 

32.6136, 37.1806, 40.2916, 

55.9914, 60.0166, 71.3698, 

89.0425, 101.2383, 


80, 1.2748, 

, 2.0876, 

5.7028, 

, 8.0464, 

12.0244, 

15.9610, 

19.6768, 20.5562, 

27.0428, 

4, 0.4778, 

, 1.4050, 

6.2150 , 

, 10.1592, 

18.4006, 

28.5824, 

40.2916, 43.4026, 


G285.60 F 
G279.20 0 


NAMELISTS 

LODE 

RKT=T , P= 36 0 . 0 , PSIA=T , OF=T , OFSKED= 3.84, 

SUPAR=30. 0,200. 0,400 .0,60 0.0, 1024. 0,ECRAT=4. 223, 

EEND 

REACTIONS 

H + OH = H20 , A=8 . 4E21 , N=2.0 , B=0., (AR) BAULCH 72 (A) 10U 

O ► H = OH , A=3 . 62E18 , N=l. , B=0., CAR) JENSEN 78 (B) 30U 

O + O = 02 , A=1.9E13 , N=0 . , B=-1.79, CAR) BAULCH 76 CA) 10U 

H + H = H2 , A=6 . 4E17 , N=l. , B=0., CAR) BAULCH 72 CA) 30U 

END TBR REAX 

H2 + OH = H20 + H , A=2.20E13, N=0.00, B=5.15, BAULCH 72 CA) 2U 

OH + OH = H20 + 0 , A=6 . 30E12 , N=0.00, B= 1.09, BAULCH 72 CA) 3U 

H2 + O = H + OH , A= 1 . 8 0E1 0 , N=-l. , B=8.9, BAULCH 72 CA) 1 . 5U 

02 + H = O + OH , A=2 . 2E1 4 , N=0. , B=16.8, BAULCH 72 CA) 1.5U 

LAST REAX 

THIRD BODY REAX RATE RATIOS 
SPECIES H2,5. ,5. ,5. ,4. , 

SPECIES H20 ,17. ,5. ,5. ,10., 

SPECIES 02,6. ,5. ,11 . ,1.5, 

SPECIES H, 12. 5, 12.5, 12. 5,25. , 

SPECIES 0,12.5,12.5,12.5,25., 

SPECIES OH, 12. 5, 12. 5, 12. 5, 25. , 

LAST CARD 
EODK 

EP= 0.0, JPRNT=-2 , MAVISP= 1 , XMC 1 ) = 1 . 0 , NJPRNT=7 , 

HI = 0 .01, HMIN=0 .01, HMAX= 0 .01, 

EEND 
ETRANS 
MP=20 0 , 

EEND 

EMOC 
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EXITPL= . FALSE . , EPU= 0 . 0 5 , 

CEND 

EBLM 

IHFLAG=0 , NTQW=1 9 , 

TQW= 126 0. 0, 126 0. 0, 126 0.0, 277 0.0, 3330.0, 3960. 0,4050.0, 

10 08.61,1039.24,886.05,610.30,571 . 14, 549. 45, 
528.76,538.35,520.14,530.12,516.42,525.62, 

XTQW= -13. 0,-6. 0,-2. 26, -2. 14, -1.57, -0.89, 0.0, 
1.8722,3.06,4.406,8.4344,13.9624, 

23.6164,32.6124,40.3644,50 .2504,62.4584,77 .9624,95.7704, 
APROF=30 . ,20 0 . ,400 .0, 1009. 0 , NPROF=4 , KDTPLT=1 , 

Kt1TPLT= 1 , KTWPLT=1 , XSEG=- 14 .0,10. 0,33. 0,56. 0,79. 0,101. 23, NSEGS=5 , 
XINOC1) =-14. 0,-12. 0,-1 0.0, -8.0, -6.0, -4.0, 

RINO(l) =2. 054, 2. 054, 2. 054, 2. 054, 2. 054, 2. 054, 

UEO(l) =8 1.0, 216. 0,351. 0,486. 0,621. 0,757.0, 

TEO(l) =5450 .0,5449.6,5448.4,5446 .2,5444.0,5442.8, 

PEO(l 1=360.0,359. 6, 358. 3, 357. 65, 356. 0,355. 9, 

NTR=700 , 

LEND 


Input for Reading 113 

LOW T CPHS 
H 

2 


ICO . 

4 968 


1 

200 . 
H2 

4 .968 

2 

2 

100 . 

6.729 


1 

200 . 
H20 

6.560 

2 

2 

100 . 

7.961 


1 

200 . 
O 

7.969 

2 

2 

100 . 

5.665 


1 

200 . 
OH 

5.433 

2 

2 

100 . 

7.798 


1 

200 . 
02 

7.356 

2 

2 

100 . 

6 . 956 


1 

200. 6.961 

END LOW T CPHS 


2 

TITLE 

DATA 

HIGH E NOZZLE STUDY: 

E=1000 AND READING 


CDATA 

ODE= 1 , CDK=1 , TDK= 0 , BLt1= 0 , IRPEAT= 0 , IOFF=6 , 

RSI= 0 . 5 , ASUB= 3 . ,1.5, NASUB = 2 , IRSTRT = 0 , 

ASUP=1 .5,2.0,30.0,200. 0 , 4 0 0 . 0 , 6 0 0 . 0 , 1 024 . 0 , NASUP=7 , 
ECRAT = 4 . 223 , RI=2 . ,THETAI = 25. ,RWTU=2. , 

ITYPE=0 ,IWALL=4,RWTD=0 . 4 , THETA=39 . 41 , 

THE = 7 . 94 , NMS=36 , 

RS= 0.0000, 1.1316, 1.1780, 1.2748, 

1.3704, 1.4198, 1.7658, 2.0876, 

2.4916, 3.0372, 5.2522, 5.7028, 

6.5606, 7.0420, 7.5072, 8.0464, 

9.0636, 9.5394, 10.6852, 12.0244, 

12.7284, 14.0056, 15.2408, 15.9610, 

16.6554, 17.3210, 18.7560, 19.6768, 20.5562, 

21.6918, 23.7554, 24.6744, 27.0428, 

29.1286, 30.1862, 32.0156, 

ZS = 0.0000, 0.3094, 0.3654, 0.4778, 

0.5908, 0.6480, 1.0430, 1.4050, 

1.8722, 2.5246, 5.5320, 6.2150, 

7.5802, 8.3862, 9.1920, 10.1592, 

12.0810, 13.0230, 15.4056, 18.4006, 

20.0684, 23.2606, 26.5588, 28.5824, 

30.6060, 32.6136, 37.1806, 40.2916, 43.4026, 

47.6334, 55.9914, 60.0166, 71.3698, 

82.7063, 89.0425, 101.2383, 

LEND 


11 


REACTANTS 

H 2. 00 100. G284.20 F 

0 2. 00 100. G277.00 0 

NAMELISTS 

CODE 

RKT = T , P = 3 56 . 9 , PSIA=T , OF=T , OFSKED=4 . 36 , 

SUP AR=30. 0,200. 0,400. 0,600. 0,1024.0, ECRAT=4 . 223 , 

LEND 

REACTIONS 

H + OH = H20 , A=8 . 4E2 1 , N=2.0 , B=0., (AR) BAULCH 72 (A) 10U 

0 + H = OH , A=3 . 62E18 , N=l. , B=0., CAR) JENSEN 78 (B) 30U 

O + O = 02 , A=1 . 9E13 , N= 0 . , B = -1.79, CAR) BAULCH 76 CA) 10U 

H + H = H2 , A = 6 . 4E1 7 , N=l. , B=0., CAR) BAULCH 72 CA) 30U 

END TBR REAX 

H2 + OH = H20 *■ H , A=2.20E13, N=0.00, B=5.15, BAULCH 72 CA) 2U 

OH + OH = H20 + O , A=6 . 30E12 , N=0.0Q, B= 1.09, BAULCH 72 CA) 3U 

H2 + 0 = H + OH , A=1 . 80E10 , N=-l. , B=8.9, BAULCH 72 CA) 1 . 5U 

02 + H = O + OH , A=2 . 2E1 4 , N=0. , B=16.8, BAULCH 72 CA) 1.5U 

LAST REAX 

THIRD BODY REAX RATE RATIOS 
SPECIES H2,5. ,5.,5. ,4. , 

SPECIES H20,17. ,5. ,5. ,10. , 

SPECIES 02,6. ,5. ,11. ,1.5, 

SPECIES H, 12. 5, 12. 5, 12. 5, 25. , 

SPECIES 0,12.5,12.5,12.5,25., 

SPECIES OH, 12. 5, 12. 5, 12. 5, 25. , 

LAST CARD 
EODK 

EP = 1 0 . 0 , JPRNT=-2 ,MAVISP=1 , XMC 1 ) =1 . 0 , NJPRNT=7 , 

HI=0 . 01 ,HNIN=0 . 01 ,HMAX=0 .01, 

EEND 
ETRANS 
MP=200 , 

EEND 

EMOC 

EXITPL=. FALSE., EPW=0. 05, 

LEND 

EBLM 

IHFLAG= 0 , NTQW=1 9 , 

TQW= 126 0.0, 126 0.0, 126 0.0, 27 7 0.0, 3330. 0,396 0.0, 40 50.0, 
1061.73,1080.00,934.61,641.91,600.62,574.11, 

5 5 1.5 0 ,^55. 18, 5 36. 6 3, 54 1.4 0,528. 52, 5 33. 44, 

XTQU1= -13. 0,-6. 0,-2. 26 , -2. 14, -1.57 , -0.89, 0.0, 
1.8722,3.06,4.406,8.4344,13.9624, 
23.6164,32.6124,40.3644,50.2504,62.4584,77 .9624,95.7704, 

APROF= 3 0 . ,20 0. ,400. 0,1009. 0 , NPROF=4 , KDTPLT=1 , NTR=7 00 , 

KMTPLT= 1 , KTWPLT= 1 , XSEG=- 1 4.0,10.0,33.0,56.0,79.0,101.23, NSEGS=5 , 
XINOC1)=-14.0, -12.0, -10.0, -8.0, -6.0, -4.0, 

RINOC 1 ) = 2 . 054, 2 . 054,2.054,2.0 54,2.054,2.054, 

UEOC 1 ) =8 0.0, 21 3. 0,346. 0,480. 0,6 13. 0,746.0, 

TEOC1) = 57 31. 8, 57 31. 0,5729. 7, 5728. 2, 57 26. 0,5724. 8, 

PEOC 11=356.9,356 .0,355.3,354.65,353.0,352.9, 

EEND 


Input for Reading 114 


LOW T 
H 

CPHS 

100 . 

4 . 968 

200 . 

4.968 

H2 


100 . 

6 . 729 

200 . 

6 560 

H20 


100 . 

7.961 

200. 

7.969 

O 


100. 

5.665 

200 . 

5.433 


2 

2 

2 

2 


1 

2 

1 

2 

1 

2 


12 


1 

2 


OH 2 

100. 7.798 1 

200. 7.356 2 

02 2 

100. 6.956 1 

200. 6.961 2 

END LOW T CPHS 

TITLE HIGH E NOZZLE STUDY: E=1000 AND READING = 114 
DATA 
CDATA 

ODE= 1 , ODK= 1 , TDK= 1 , BLM= 1 , IRPEAT=2 , IOFF=6 , 

RSI=0 . 5 , ASUB=3 .,1.5, NASUB = 2 , IRSTRT=0 , 

ASUP= 1.5*2.0,30. 0,200. 0,400. 0,600. 0,1024.0, NASUP=7 , 
ECRAT=4 . 22 3 , RI=2 . ,THETAI=25. ,RWTU=2. , 

ITYPE-0 , IWALL=4 , RWTD=0 . 4 , THETA=39 .41# 

THE=7 . 94 , NWS =3 6 , 

RS= 0.0000, 1.1316, 1.1780, 1.2748, 

1.3704, 1.4198, 1.7658, 2.0876, 

2.4916, 3.0372, 5.2522, 5.7028, 

6.5606, 7.0420, 7.5072, 8.0464, 

9.0636, 9.5394, 10.6852, 12.0244, 

12.7284, 14.0056, 15.2408, 15.9610, 

16.6554, 17.3210, 18.7560, 19.6768, 20.5562, 

21.6918, 23.7554, 24.6744, 27.0428, 

29.1286, 30.1862, 32.0156, 

ZS= 0.0000, 0.3094, 0.3654, 0.4778, 

0.5908, 0.6480, 1.0430, 1.4050, 

1.8722, 2.5246, 5.5320, 6.2150, 

7.5802, 8.3862, 9.1920, 10.1592, 

12.0810, 13.0230, 15.4056, 18.4006, 

20.0684, 23.2606, 26.5588, 28.5824, 

30.6060, 32.6136, 37.1806, 40.2916, 43.4026, 

47.6334, 55.9914, 60.0166, 71.3698, 

82.7063, 89.0425, 101.2383, 


80, 1.2748, 

I, 2.0876, 

5.7028, 

!, 8.0464, 

12.0244, 

15.9610, 

19.6768, 20.5562, 

27.0428, 

>4, 0.4778, 

), 1.4050, 

6.2150, 

), 10.1592, 

, 18.4006, 

, 28.5824, 


G283.90 F 
G275.80 O 


NAMELISTS 

LODE 

RKT=T,P=360 . 9, PSIA=T , OF=T , OFSKED=5 . 08 , 

SUPAR=30 . 0 ,200 . 0,400 .0,600 . 0, 1024. 0,ECRAT=4 .223, 

SEND 

REACTIONS 

H + OH = H20 , A=8 . 4E21 , N=2.0 , B=0., (AR) BAULCH 72 (A) 10U 

O + H = OH . A=3 . 62E18 , N=l. , B=0., CAR) JENSEN 78 CB) 30U 

O + O = 02 , A= 1 . 9E1 3 , N=Q . , B=-1.79, (AR) BAULCH 76 (A) 10U 

H + H = H2 , A=6 . 4E17 , N=l. , B=0., (AR) BAULCH 72 (A) 30U 

END TBR REAX 

H2 + OH = H20 + H , A=2.20E13, N=0.00, B=5.15, BAULCH 72 (A) 2U 

OH + OH = H20 + O , A=6 . 30E12 , N=0.00, B= 1.09, BAULCH 72 (A) 3U 

H2 + O = H + OH , A=1 .80E10 , N=-l. , B=8.9, BAULCH 72 (A) 1 . 5U 

02 + H = O + OH , A-2 . 2E14 , N=0. , B=16.8, BAULCH 72 (A) 1 . 5U 

LAST REAX 

THIRD BODY REAX RATE RATIOS 
SPECIES H2,5. ,5. ,5. ,4. , 

SPECIES H20, 17 . ,5. ,5. , 10 . , 

SPECIES 02,6 . , 5 . , 11 . , 1 . 5, 

SPECIES H, 12. 5, 12. 5, 12. 5, 25. , 

SPECIES 0,12.5,12.5,12.5,25., 

SPECIES OH, 12. 5, 12. 5, 12. 5, 25. , 

LAST CARD 
£ODK 

EP=0 . 0 , JPRNT=-2 ,MAVISP=1 ,XM( 1 )=1 . 0 ,NJPRNT=7 , 

HI=0 . 01 ,HMIN=0 .01 ,HMAX=0 .01, 

CEND 
ETRANS 
MP=20 0 , 

CEND 

CMOC 

EXITPL=. FALSE. ,EPW=0. 05, 

CEND 
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Or FCiuK QUALITY 


CBLM 

IHFLAG=0,NTQW=17, 

TQW= 1 26 0 . 0,2770.0,3330.0,396 0.0,4050.0, 

1169. 58, 11 28. 31, 972. 62, 663. 76, 618. 64, 58 9. 39, 
567.31,569.98,548.78,549.66,535.97,537.46, 

XTQW= -2.26,-2.14,-1.57,-0.89,0.0, 
1.8722,3.06,4.406,8.4344,13.9624, 
23.6164,32.6124,40.3644,50.2504,62.4584,77.9624,95.7704, 
APROF=30 . ,20 0 .,400.0,1009.0, NPROF=4 ,KDTPLT=1 , 

KMTPLT=1 , KTWPLT= 1 , XSEG=- 1 4.0,10.0,33.0,56.0,79.0,101.23, NSEGS=5 , 
XINOC 1)=-14.0, -12.0, -10.0, -8.0, -6.0, -4.0, 


RINO(l 1=2.054,2.054, 2. 054, 2. 054, 2. 054, 2. 054, 

UEO( 1) = 35. 0,165. 0,300. 0,425. 0,56 0.0, 685.0, 

TEO( 1 1=6014.8,6013. 6, 6012. 4, 6011 .2,6010 . 0,6008.8, 
PEO( 11=360. 7, 360 .0,359.3,358.65,358.0,357.3, 
NTR=80 0 , 

GEND 


Input for Reading 115 


LOW T 
H 

CPHS 

2 


100 . 

4.968 


1 

200 . 
H2 

4 . 968 

2 

2 

100 . 

6 .729 


1 

200 . 
H20 

6.560 

2 

2 

100 . 

7.961 


1 

200 . 
O 

7.969 

2 

2 

100 . 

5.665 


1 

200 . 
OH 

5.433 

2 

2 

100 . 

7.798 


1 

200 . 
02 

7 . 356 

2 

2 

100. 

6 . 956 


1 

200. u.961 

END LOW T CPHS 


2 

TITLE 

DATA 

HIGH E NOZZLE 

: STUDY: 

E=1000 AND READING 


GDATA 

ODE= 1 , ODK= 1 , TDK=1 , BLM=1 , IRPEAT=2 , IOFF=6 , 
RSI=0 . 5. ASUB=3. , 1 . 5 , NASUB=2 , IRSTRT=2 , 


ASUP=1 .5,2.0,30. 0,200.0,400.0,600.0, 1 024 . 0 , NASUP=7 , 
ECRAT=4 . 223 , P.I=2 . , THETAI=25 . ,P.WTU=2. , 

ITYPE=0 , IWALL=4 , RWTD=0 . 4 ,THETA=39 . 41 , 

THE = 7 . 94 , NWS=36 , 

RS = 0.0000, 1.1316, 1.1780 

1.3704, 1.4198, 1.7658, 

2.4916, 3.0372, 5.2522, 

7.5072, 

10 . 6852 , 

15.2408, 

18.7560, 19.6768, 20.5562, 


7.0420, 
9.5394, 
14.0056 , 
17.3210, 


6.5606 
9.0636 , 

12.7284, 

16.6554, 

21 .6918, 

29.1286, 

ZS= 0 . 0000 , 

0.5908, 0.6480 

1.8722, 2.5246, 

7.5802, 8.3862 


12.0810, 
20 .0684, 
30.6060, 
47.6334, 
82.7063, 
SEND 

REACTANTS 
H 2 . 

0 2 . 


23.7554, 24.6744, 

30.1862, 32.0156, 

0.3094, 0.3654, 

1.0430, 

5 . 5320 , 
9.1920, 
13.0230, 15.4056, 

23.2606, 26.5588, 

32.6136, 37.1806, 40 

55.9914, 60.0166, 71 

89.0425, 101.2383, 


1.2748 
2. 0876, 
5.7028, 
8.0464, 
12.0244, 
15.9610, 
19.6768, 
27.0428, 


0.4778, 
1.4050, 
6.2150, 
10 .1592, 
18.4006, 
28.5824, 
.2916 
. 3698 


43.4026, 


00 

00 


100 . 
100 . 


G283.50 F 
G275.90 0 
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NAMELISTS 

CODE 


RKT=T,P=355. 3,PSIA=T,OF=T,OFSKED=5.49, 

SUPAR=3Q. 0,200. 0,400. 0,600. 0,102 4.0, ECRAT=4 .223, 

SEND 

REACTIONS 

H + OH = H20 , A = 8 . 4E21 , 

0 + H - OH , A=3 . 62E18 

0 + O = 02 , A=1 . 9E1 3 , 

H + H = H2 , A = 6 . 4E17 , 


END TBR REAX 

H2 + OH = H20 + H , 
OH + OH = H20 + O 

H2 + O = H + OH 
02 + H = O + OH 
LAST REAX 


A=2 . 20E1 3 , 

, A = 6 . 3 0E1 2 , 
A=1 .80E10, 
A=2 . 2E1 4 , 


N=2 . 0 , B=0 . , CAR) BAULCH 72 (A) 10U 

, N=l. , B=G . , CAR) JENSEN 78 (B) 30U 
N=0. , B=-1.79, CAR) BAULCH 76 CA) 10U 

N=l. , B=0., CAR) BAULCH 72 CA) 30U 

N = 0 . 0 0 , B=5 . 15 , BAULCH 72 CA) 2U 

N= 0 . 0 0 , B= 1.0 9, BAULCH 7 2 CA) 3U 
N=-l. , B=8.9, BAULCH 72 (A) 1 . 5U 

N=0. , B=16 . 8 , BAULCH 72 (A) 1.5U 


THIRD BODY REAX RATE RATIOS 
SPECIES H2,5. ,5. ,5. ,4. , 

SPECIES H20, 1 7 . ,5. ,5. ,10. , 

SPECIES 02,6. ,5. ,11. ,1.5, 

SPECIES H, 12. 5, 12. 5, 12. 5, 25. , 

SPECIES 0,12.5,12.5,12.5,25., 

SPECIES OH, 12.5, 12.5,12.5,25. , 

LAST CARD 
BOOK 

ET=0 . 0 , JPRNT=-2,MAVISP=1 ,XMC1)=1 . 0,NJPRNT=7, 

HI = 0 .01, HMIN= 0 .01, HMAX=0 .01, 

BEND 

STRANS 

MP = 2 0 0 , 

BEND 

BMOC 

EXITPL=. FALSE. ,EPW=0. 05, 

BEND 

BBLM 

IHFLAG=0 ,NTQW=19, 

TQW =1260. 0,1260. 0,1260. 0,2770. 0,3330. 0,3960.0, 4050.0, 
1156.32,1117.60,985.76,682.10,637.87,603.26, 

58 1.66, 57 9. 76, 559. 11, 557. 7 1,542. 84, 542. 02, 

XTQW = -13. 0,-6. 0,-2. 26, -2. 14, -1.57, -0.89, 0.0, 

1.8722,3.06,4.406,8.4344,13.9624, 
23.6164,32.6124,40.3644,50.2504,62.4584,77.9624,95.7704, 

APROF= 3 0 . , 200. ,400.0,1009.0, NFROF = 4 , KDTPLT=1 , NTR=800 , 

KMTPLT= 1 , KTWPLT= 1 , XSEG=- 14 . 0, 10. 0,33. 0,56.0, 79. 0,I0t. 23, NSEGS=5 , 
XINOC1) =-14. 0,-12. 0,-10. 0,-8. 0,-6. 0,-4 .0, 

RINOC 1) = 2. 0 54, 2. 054, 2. 0 54, 2. 0 54, 2. 054, 2. 054, 

UEOC 1 ) =8 0 . 0,213. 0,346. 0,48 0.0, 613. 0,746.0, 

TEOCl)=6 124 .7,6123.8,6122.0,6120.2,6119.0,6118.0, 

PEOC1) =355. 3, 355. 0,354. 3, 353. 65, 352. 9, 352.0, 

BEND 


Input for Reading 120 


LOW 

T CPUS 


H 


2 

100 . 

4 . 968 


200 . 

4.968 


M2 


2 

100 , 

6.729 


200 , 

6.560 


H20 


2 

100 . 

7.961 


200 . 

7.969 


0 


2 

100 . 

5.665 


200 . 

5.433 


OH 


2 

100 . 

7.798 


200 . 

7 . 356 


02 


2 


1 

2 

1 

2 

1 

2 

1 

2 

1 

2 


15 



1 

2 


100. 6.956 

200. 6.961 

END LOW T CPHS 
TITLE HIGH E NOZZLE STUDY: E=1000 AND READING = 120 
DATA 
EDATA 

OD£=l ,ODK=l ,TDK=0,BLM=0 ,IRPEAT=0 ,IOFF = 6 , 

RSI = 0 . 5 , ASUB=3 . , 1 . 5 ,NASUB = 2 ,IP.STRT=0 , 

ASUP=1 .5, 2. 0,30. 0,200. 0,400. 0,600. 0.1029.0, NASUP=7 , 

ECP.AT = 4 . 22 ^ . RI=2 . ,THETAI=25. ,RWTU-2. , 

ITYPE=0 ,IWALL=4 ,RWTD=0 .4 ,THETA=39 . 41 , 

THE=7 . 94,NWS=36 , 

RS= 0.0000, 1.1316, 1.1780, 1.2748, 

1.3704, 1.4198, 1.7658, 2.0876, 

2.4916, 3.0372, 5.2522, 5.7028, 

6.5606, 7.0420, 7.5072, 8.0464, 

9.0636, 9.5394, 10.6852, 12.0244, 

12.7284, 14.0056, 15.2408, 15.9610, 

16.6554, 17.3210, 18.7560, 19.6768, 20.5562, 

21.6918, 23.7554, 24.6744, 27.0428, 

29.1286, 30.1862, 32.0156, 

ZS= 0.0000, 0.3094, 0.3654, 0.4778, 

0.5908, 0.6480, 1.0430, 1.4050, 

1.8722, 2.5246, 5.5 32 0 , 6.2150, 

7.5802, 8.3862, 9.1920, 10.1592, 

12.0810, 13.0230, 15.4056, 18.4006, 

20.0684, 23.2606, 26.5588, 28.5824, 

30.6060, 32.6136, 37.1806, 40.2916, 43.4026, 

47.6334, 55.9914, 60.0166, 71.3698, 

82.7063, 89.0425, 101.2383, 

EEND 

REACTANTS 

H 2. 00 100. G294.40 F 

0 2. 00 100. G287.50 0 

NAMELISTS 

CODE 

RKT = T,P = 355 . 2 , PSIA = T , OF = T , OFSKED=4 . 30 , 

SUPAR=30. 0,200. 0,400. 0,600. 0,1024.0, ECRAT=4 .223, 

EEND 

REACTIONS 

H + OH = H20 , A=8 . 4E2 1 , N=2.0 , B=0., (AR) BAULCH 72 (A) 10U 

O + H = OH , A=3 . 62E18 , N=l. , B=0., (AR) JENSEN 78 (B) 30U 

O + O = 02 , A=1 . 9E1 3 , N=0. , B=-1.79, (AR) BAULCH 76 (A) 10U 

H + H = H 2 , A=6 . 4E1 7 , N=l. , B = 0., (AR) BAULCH 72 (A) 30U 

END TBR REAX 

H2 + OH = H20 + H , A=2.20E13, N=0.00, B=5.15, BAULCH 72 (A) 2U 

OH J . OH = H20 + O , A = 6 . 30E12 , N=0.00, B= 1.09, BAULCH 72 (A) 3U 

H2 + O = H + OH , A= 1 . 80E1 0 , N=-l. , B=8.9, BAULCH 72 (A) 1 . 5U 

02 + H = 0 + OH , A=2 . 2E1 4 , N=0. , B=16.8, BAULCH 72 (A) 1 . 5U 

LAST REAX 

THIRD BODY REAX RATE RATIOS 
SPECIES H2,5. ,5. ,5. ,4. , 

SPECIES H20, 17. ,5. ,5. , 10 . , 

SPECIES 02,6 . ,5. ,11. ,1 .5, 

SPECIES H, 12. 5, 12. 5, 12. 5, 25. , 

SPECIES 0,12.5,12.5,12.5,25., 

SPECIES OH, 12. 5, 12. 5, 12. 5, 25. , 

LAST CARD 
CODE 

EP=10 . 0 , JPRNT=-2 ,MAVISP=1 ,XM( 1 )=1 . 0 ,NJPRNT=7 , 

HI=0 . 01 ,HMIM=0 . 01 ,HMAX=0 . 01 , 

EEND 
ETRANS 
MF = 2 0 0 , 

LEND 

CMOC 

EXITPL= . FALSE . , EPW= 0 .05, 

LEND 

CBLM 

IHFLAG= 0 , NTQW=1 9 , 

TQw= 126 0.0 , 126 0.0 , 126 0.0, 27 7 0.0, 3330. 0,396 0.0, 4050.0 , 
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1149. 34, 9 58. 00, 863. 83, 635. 66, 548. 76, 577. 38, 
565.98,559.23,551.66,550.73,543.98,538.47, 

XTQW= - 13 . 0,-6. 0,-2. 26, -2. 14, -1.57, -0.89, 0.0, 
1.8722,3.06,4.406,8.4344,13.9624, 

23 6164 , 32 . 6124 , 40. 36 44 , 50. 25 04, 62 . 4 584, 77. 96 24 , 95. 7704, 
APROF-30 . ,200 . ,400.0,1009.0, NPROF=4 , KDTPLT=1 , NTR=800, 

KMTPLT= 1 ,KTWPLT=1 , XSEG=- 14 . 0,10. 0,33. 0,56. 0,79. 0,101. 23, NSEGS-5, 
XIHOC 1)=- 14.0, -12.0, -10.0, -8.0, -6.0, -4.0, 

RIMOC 1 5=2.054, 2. 054, 2. 054, 2. 054, 2. 054, 2. 054, 

UEO( 1 )=8Q . 0,213. 0,346. 0,480. 0,613. 0,7 46.0, 

TEO (15 = 5709. 6, 5709. 0,5707. 5, 5706. 0,5704. 0,5702. 2, 

PECK 15 = 355. 2, 355. 0,354. 3, 353. 65, 3 52. 9, 352.0, 

£,END 


Input for Reading 121 


LOW T CPHS 

H 2 

100. 4.968 1 

200. 4.968 2 

H2 2 

100. 6.729 1 

200. 6.560 2 

H20 2 

100. 7.961 1 

200. 7.969 2 

0 2 

100. 5.665 1 

200. 5.433 2 

OH 2 

100. 7.798 1 

200. 7.356 2 

02 2 

100. 6.956 1 

200. 6.961 2 

END LOU T CPHS 

TITLE HIGH E NOZZLE STUDY: E=1000 AND READING = 121 
DATA 
SDATA 

ODE=l , ODK=l , TDK=1 , BLM=1 , IRPEAT=2 , IOFF=6 , 

RSI= 0.5, ASUB = 3 . ,1.5, NASUB=2 , IRSTRT=0 , 

A3UP=1 . 5,2.0,30.0,200.0,400.0,600.0,1024.0 ,NASUP=7 , 
ECRAT=4 . 223 , RI=2 . ,THETAI=25. ,RWTU=2. , 

ITYPE= 0 , IWALL=4 , RWTD=0 . 4 ,THETA = 39 . 41 , 


THE=7 . 94 ,NWS=36 , 

RS= 0.0000, 1.1316, 1.1780, 1.2748, 

1.3704, 1.4198, 1.7658, 2.0876, 


2.4916 , 

3.0372, 5.2522, 

5 .7028, 

6.5606 

, 7.0420, 7.5072, 

8 .0464, 

9.0636, 

9.5394, 10.6852, 

12.0244, 

12.7284, 

14.0056, 15.2408, 

15.9610, 

16.6554, 

17.3210, 18.7560, 19 

.6768, 20.5562 

21.6918, 

23.7554, 24.6744, 27.0428, 

29.1286 , 

30.1862, 32.0156, 


ZS= 0.0000 

, 0.3094, 0.3654, 

0.4778, 

0.5908 

, 0.6480, 1.0430, 

1.4050, 

1.8722, 

2.5246, 5.5320, 

6.2150, 

7.5802 

, 8.3862, 9.1920, 

10.1592, 

12.0810, 

13.0230, 15.4056, 

18.4006, 

20.0684, 

23.2606, 26.5588, 

28.5824, 

30.6060, 

32.6136, 37.1806, 40 

.2916, 43.4026 

47.6334, 

55.9914, 60.0166, 71 

.3698, 

82.7063, 

89.0425, 101.2383, 


EEND 



REACTANTS 



H 2. 


00 100. 

O 2. 


00 100. 


G295.00 F 
G288.30 O 
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NAMELISTS 

CODE 


RKT=T,P=360 . 0,PSIA=T,OF=T,OFSKED=5 . 1 1 
SUPAR= 30. 0,20 0.0, 40 0.0, 6 00. 0,1 024.0, 
SEND 


REACTIONS 




H + 

OH = 

H20 

, A=8 . 4E2 1 , 

, N=2 . 0 

O + 

H = 

OH 

, A=3 . 62E18 

, N = l. 

O + 

O = 

02 

, A=1 . 9E1 3 , 

, N = 0. 

H + 

H = 

H2 

, A=6 . 4E1 7 , 

, N= 1 . 

END TBR REAX 



H2 

+ OH 

= H20 + H , 

A=2 . 2 0E1 3 , 

N= 0 . 0 

OH 

+ OH 

= H20 + O 

, A=6 . 30E12 , 

, N = 0 . 

H2 

+ O = 

H + OH 

A= 1 . 8 0E1 0 , 

N = -l . 

02 

+ H = 

O + OH 

A=2 . 2E14 , 

N= 0 . 


LAST P.EAX 

THIP.D BODY REAX RATE RATIOS 
SPECIES H2,5. ,5. ,5. ,4. , 

SPECIES H20 , 17 . ,5. ,5 . , 10 . , 
SPECIES 02,6. ,5. ,11. ,1 .5, 
SPECIES H, 12. 5, 12. 5, 12. 5, 25. , 
SPECIES 0,12.5,12.5,12.5,25., 
SPECIES OH, 12. 5, 12. 5, 12. 5, 25., 
LAST CARD 
EODK 


ECRAT=4 . 223 , 


, B=0 . , (AR) BAULCH 72 (A) 10U 
, B=0 . , (AR) JENSEN 78 (B) 30U 
, B=-1.79, (AR) BAULCH 76 (A) 10U 
, B=0 . , (AR) BAULCH 72 (A) 30U 

0, B=5.15, BAULCH 72 (A) 2U 

00, B= 1.09, BAULCH 72 (A) 3U 
. B=8 . 9 , BAULCH 72 (A) 1.5U 

, B= 16 . 8 , BAULCH 72 (A) 1.5U 


EP= 0 . 0 , JPP.NT=-2 ,MAVISP=1 ,XM( 1 ) = 1 . 0 ,NJPRNT = 7 , 


HI = 0 . 0 1 , HMIN= 0 .01, HMAX = 0 .01, 

CEND 
ETRANS 
MP = 20 0 , 

LEND 

EMOC 

EXITPL=. FALSE. ,EPW=0 .05, 

LEND 

EBLM 

IHFLAG=0 , NTQU= 1 9 , 

TQW= 126 0.0, 126 0.0, 126 0.0, 27 7 0.0, 3 33 0.0, 3 96 0.0, 4050.0, 
1057.56,922.32,873.83,641.26,598.97,572.43, 
563.59,558.51,552.33,545.82,542.90,537.82, 

XTQW= -13. 0,-6. 0,-2. 26, -2. 14, -1.57, -0.89, 0.0, 
1.8722,3.06,4.406,8.4344,13.9624, 
23.6164,32.6124,40.3644,50.2504,62.4584,77.9624,95.7704, 

APROF= 3 0 . , 200. ,400.0,1009.0, NPROF=4 ,KDTPLT=1 , 

KMTPLT= 1 , KTWPLT = 1 , XSEG= - 14.0,10.0,33.0,56.0,79.0,101.23, NSEGS= 5 , 
XINO(l)=- 14. 0,-12. 0,-1 0.0, -8.0, -6.0, -4.0, 

RIMO(l) =2. 054, 2. 054, 2. 054, 2. 054, 2. 054, 2. 054, 

UEO(l) =8 1.0,21 6. 0,3 5 1.0, 48 6. 0,621. 0,757.0 , 

TEO(l) = 56 14. 0,56 12. 6, 56 10. 4, 56 08. 2, 56 07. 0,56 05. 8, 

PEO(l) = 36 0.0, 359. 6, 358. 3, 357. 65, 356. 0,355. 9, 

NTR = 80 0 , 

SEND 

/EOF 


Input for Reading 137 


LOW T CPHS 
H 

100 . 

200 . 

H2 
100 . 

200 . 

H20 
100 . 

200 . 

O 

100 . 

200 . 

OH 
100 . 


4 . 968 
4 . 968 

2 

6 . 729 
6 . 560 

2 

7 . 961 
7 . 969 

2 

5.665 

5.433 

2 

7 .798 


1 

2 

1 

2 

1 

2 

1 

2 

1 
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200. 7.356 2 

02 2 

100. 6.956 1 

200. 6.961 2 

END LOW T CPHS 

TITLE HIGH E NOZZLE STUDY: E=1000 TRUNCATED TO 900:1 AND READING = 137 
DATA 
CDATA 

ODE=l ,ODK=l,TDK=l,BLM=l ,IRPEAT=2, 

RSI=0.5, ASUB = 3. , 1 . 5 , NASUB=2 , IOFF=6 , IRSTRT=2 , 

ASUP=1 . 5,2.0.30.0,200.0,930. 0 ,NASUP=5 , 

ECP.AT-9 . 22 3 , RI = 2 . ,THETAI=25. ,RWTU=2. , 

ITYPE=0 ,IWALL=9 .RWTD-0 . 9 ,THETA=39 . 91 , 


THE=15 . 5 ,NWS=28 , 


RS= 0.0000, 
1.2986, 
2.7939, 

6 . 0729, 
9.3032, 
13.7082, 
17.6959, 

ZS = 0.0000 

0.5058, 

2. 1699, 
6.7929, 
12.5520 , 
22.9989, 
33.6179, 
EEND 

REACTANTS 

H 2. 

0 2 . 


1.1109, 1 

1.3962, 1. 

3.3786, 9. 

6.8039, 7. 

9.7722, 11. 

19.2982, 19. 

19.2220, 20. 

, 0.2819, 

0.5620, 0. 

2.9508, 9. 

7.9832, 8. 

13.9938, 16. 

29.0228, 25. 

38.7360, 91. 


1 

1 

5 

7 

11 

16 

20 


.1592, 
9996 , 
7169, 

2766 , 
0168, 

8709 , 

1216 , 

0 . 3379, 
6766 , 0 

7509 , 
7892, 

1258 , 
5972, 

8970 , 


6 

9 

17 

29 

99 


.2510, 

.5226, 

.5916, 

.7792, 

.6616, 

.3112, 

.8512, 

0 . 9996 
.7662, 
.0992, 
.6756 , 

. 5666 , 
.5992, 
.9869, 


00 

00 


100 . 
100 . 


G291 .20 
G285 .00 


F 

0 


NAMELISTS 

CODE 

PKT=T,P=356 .8, PSIA=T , OF=T, OFSKED=9 . 29 , 
SUPAR=30 .0,200.0,930. 0 ,ECRAT=9 .223, 
EEND 


REACTIONS 


H + OH = H20 

, A=8 . 9E2 1 

O + H = OH 

, A= 3 . 62E18 

O+O =02 

, A= 1 . 9E1 3 

H + H = H2 

, A=6 . 9E1 7 


, B=0., (AR) BAULCH 72 (A) 10U 
, B-0 . , (AR) JENSEN 78 (B) 30U 
, B=-l . 79 , CAR) BAULCH 76 (A) 10U 
, B=0., (AR) BAULCH 72 (A) 30U 


H2 + OH = H20 + H , A=2.20E13, N=0.00, B=5.15, BAULCH 72 (A) 2U 

OH + OH = H20 + O , A=6.30E12, N-0.00, B = 1.09, BAULCH 72 (A) 3U 

H2 + 0 = H + OH , A= 1 . 8 0E1 0 , N=-l. , B = 8.9, BAULCH 72 (A) 1 . 5U 

02 + H = 0 + OH , A = 2 . 2E19 , N=0. , B=16.8, BAULCH 72 (A) 1.5U 

LAST REAX 

THIRD BODY REAX RATE RATIOS 
SPECIES H2,5. ,5. ,5. ,9. , 

SPECIES H20,17. ,5. ,5. ,10. , 

SPECIES 02,6. ,5. ,11. ,1.5, 

SPECIES H, 12. 5, 12. 5, 12 . 5,25 . , 

SPECIES 0,12.5,12.5,12.5,25., 

SPECIES OH, 12. 5, 12. 5, 12. 5, 25., 


LAST CARD 
CODE 

EP= 10 . 0 ,MAVISP=1 ,XM( 1 )=1 . 0 , JPRNT=-2, 

HI=0 . 010, HMIN= 0.010, HMAX= 0 . 010 , ARPRNT ( 1 ) = 939 . 0,NJPRNT=1, 
EEND 
ETRANS 

MP = 200 , 

EEND 

EMOC 

EXIT PL= . FALSE . , EPW= 0 .05, 

EEND 

EBLM 

IHFLAG=0 , NTQW=15 , 

TQW= 1260.0,1260.0,1260.0,2700.0,3330.0,3960.0,9050.0, 
1199.39, 958.00,863.83,635.66,598.76,577.38,565.98, 
559.23, 
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XTQW = -13. 0,-6. 0,-2. 26, -2. 14,-1 .57 ,-0.89,0.0 , 
1.8722,3.06,4.406,8.4344,13.9624, 

23.6164,32.6124,40.3644, 

APROF= 3 0 . ,200. ,396.0, HPROF= 3 , KDTPLT= 1 , NTR=800 , 

Kt1TPLT= 1 , KTWPLT= 1 , XSEG= - 14.0,-2.0,10.0,22.0,34.0,44.2, NSEGS=5 , 
XINO = - 14.0, -12.0, -10.0, -8.0, -6.0, -4.0, 

RINO = 2. 054, 2. 054, 2. 054, 2. 054, 2. 0 54, 2. 0 54, 

UEO= 85.0,227.0,368.0,510.0,652.0,793.0, 

PEO= 356.4,355.7,355.0,354.3,353.6,352.9, 

TEO= 5702.0,5699.6,5698.5,5697.2,5695.9,5695.6, 

CEND 

Appendix B 

Pressure Integration Calculation 


Pressure Integration Procedure 

The force acting on the surface of a rocket nozzle can be 
represented as the result of normal and tangential forces. The 
normal force per unit area is defined as pressure. Conversely, 
tangential force per unit area is termed shear stress. For 
symmetric nozzles, only the axially directed force is of concern 
with regard to thrust. 

According to boundary -layer theory, the normal forces 
acting on a rocket nozzle are independent of tangential forces 
(except for minor boundary-layer displacement corrections). 
As such, it is possible to experimentally distinguish inviscid 
(normal) forces from shear forces. The inviscid thrust between 
discrete nozzle area ratios can be determined by integrating 
pressure with respect to the normal surface area: 

r ^ 2 

A Thrust = I P w dA N (1) 


between the inviscid and net thrust (or thrust coefficient) is 
then the shear, or drag, decrement. 

The integration of equation (2) was carried out by performing 
a piecewise integration of measured pressures. The relationship 
between pressure and area ratio was assumed to have the form 


This form is considered accurate due to the nearly linear nature 
of the pressure versus area ratio plot when shown on log-log 
scales (figs. 9 and 10). Coefficients a and b are determined 
from measurements of pressure at two distinct area ratios. 
Usually these points represent the two end points in the 
piecewise integration. The exception to this rule is the inter- 
polation done at the lowest area ratios and the extrapolation 
performed out to the exit area ratio. 

Determination of Inviscid Thrust Gain 


where 

P w static wall pressure 
A n normal surface area 

or, in terms of thrust coefficient and area ratio, 


Acy v — 



de 


(2) 


where 

C F V vacuum thrust coefficient 
e area ratio 

P c e effective chamber pressure 

Furthermore, the net gain in thrust (or thrust coefficient) 
between area ratios can be determined by testing a nozzle 
contour truncated at various area ratios. The difference 


As previously discussed, the thrust gain of a rocket engine 
between two given area ratios is the net result of the inviscid 
thrust gain and the drag decrement. Only the axial thrust is 
of concern as radial components of thrust are assumed to 
cancel. To determine the force of the rocket engine as the result 
of normal stresses (pressure), one considers only the normal 
component of the nozzle surface area, or 

r^2 

A Thrust = I P w dA N (1) 

•W, 

and, since throat area A, is constant, 


dAft — At d€ 


Also, by definition 


ACy v — 


A Thrust 

Pc,e A t 
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Therefore, 



As discussed previously, the plot of pressure ratio versus 
area ratio is nearly linear when presented on log-log scales. 
Therefore, since equations of the form y = a(x b ) appear as 
straight lines on log-log paper, this relation was assumed; that 
is, 



Piecewise integration of equation (2) is then 


ACp |/ — 


ae 


fc+i 


b+ 1 


Table IV shows a sample calculation of the inviscid thrust 
gain. 


TABLE IV. -SAMPLE CALCULATION 
[Nozzle contour, 1000:1; area ratio, 1030; mixture ratio, 3.84.] 


Area ratio 

Pressure 

ratio, 

PJPc, 

Integration region 

Parameter 

Vacuum 

thrust 

coefficient, 

Cp.v 



a 

b 

388.01 

0.0002294 

427.50 to 499.97 

0.54749 

-1.304734 

0.01186 

499.97 

.0001648 

499.97 to 635.04 

.47666 

-1.282441 

.01905 

635.04 

.0001213 

635.04 to 799.98 

1.7719 

-1.485888 

.01682 

799.98 

.0000861 

799.98 to 974.94 

.02019 

-.816495 

.01387 

974.94 

.0000732 

974.94 to 1025.0 

“02019 

a -. 816495 
Total 

.00359 

.06519 


indicates extrapolated value. 


Integration Uncertainty 


The constants a and b are determined by using two data points 
of pressure and area ratio, where 



and 


The uncertainty in performing piecewise integration was 
examined. This was done by using the TDK program. The 
pressure output, at area ratios corresponding (approximately) 
to those at which the experimental pressure measurements were 
taken, was tabulated. The incremental thrust coefficient gain 
from area ratios of 427.5 to 1030 was predicted at 0.0600 by 
using the integration procedure described previously. This 
compares with the computationally predicted gain of 0.0610. 
Therefore, it seems reasonable to conclude that the accuracy 
of the integration procedure is within 2 percent of the actual 
conditions. 
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